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Abstract 

We measure the luminosity and color dependence and the redshift evolution of galaxy clustering 
_ in the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Ninth Data Release. 

rS . We focus on the projected two-point correlation function (2PCF) of subsets of its CMASS sample, 

■ which includes about 260,000 galaxies over ~3,300deg^ in the redshift range 0.43 < z < 0.7. To 

minimize the selection effect on galaxy clustering, we construct well-defined luminosity and color 
subsamples by carefully accounting for the CMASS galaxy selection cuts. The 2PCF of the whole 
CMASS sample, if approximated by a power-law, has a correlation length of ro = 7.93 ± 0.06/i^^Mpc 
and an index of 7 = 1.85 ± 0.01. Clear dependences on galaxy luminosity and color are found for 
the projected 2PCF in all redshift bins, with more luminous and redder galaxies generally exhibiting 
stronger clustering and steeper 2PCF. The color dependence is also clearly seen for galaxies within 
the red sequence, consistent with the behavior of SDSS-II main sample galaxies at lower redshifts. 
At a given luminosity (k-l-e corrected), no significant evolution of the projected 2PCFs with redshift 
is detected for red sequence galaxies. We also construct galaxy samples of fixed number density at 
different redshifts, using redshift-dependent magnitude thresholds. The clustering of these galaxies in 
the CMASS redshift range is found to be consistent with that predicted by passive evolution. Our 
measurements of the luminosity and color dependence and redshift evolution of galaxy clustering will 
allow for detailed modeling of the relation between galaxies and dark matter halos and new constraints 
on galaxy formation and evolution. 

Subject headings: cosmology: observations — cosmology: theory — galaxies: distances and redshifts 
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— galaxies: halos — galaxies: statisties — large-seale strueture of Universe 



1. INTRODUCTION 

Galaxies in the universe are observed to display a wide 
range of properties, such as luminosity, color, stellar 
mass, age, morphology, and spectral type. These prop- 
erties encode information about galaxy formation and 
evolution and are related to the environment hosting the 
galaxies. Different populations of galaxies are thus ex- 
pected to trace the underlying dark matter distribution 
in different ways. 

Contemporary galaxy redshift surveys, most notably 
the Sloan Digital Sky Survey (SDSS; York et al. 2000), 
have transformed the study of large-scale structure, 
enabling pristine measurements and detailed studies of 
galaxy clustering. Galaxy clustering provides a pow- 
erful approach to probe the complex relation between 
galaxies and the underlying dark matter distribution 
(e.g.. Kaiser 1984). The dependence of galaxy clustering 
on galaxy properties has been observed in numerous 
galaxy surveys (e.g., Davis & Geller 1976; Davis et al. 
1988; Hamihon 1988; Lovcday et al. 1995; Benoist et al. 
1996; Guzzo ct al. 1997; Norberg et al. 2001, 2002; 
Zehavi et al. 2002, 2005b, 2011; Budavari ct al. 
2003; Madgwick et al. 2003; Li et al. 2006; Coil et al. 
2006, 2008; Meneux et al. 2006, 2008; Wake ct al. 
2008; Swanson ct al. 2008; Meneux et al. 2009; 
Ross & Brunner 2009; Skibba et al. 2009; Loh et al. 
2010; Ross et al. 2010, 2011a; Wake et al. 2011; 
Christodoulou et al. 2012; Mostek et al. 2012). 

In general, more luminous and redder galaxies are 
found to be more strongly clustered than their fainter 
and bluer counterparts. Similarly, early-type (ellip- 
tical) galaxies exhibit stronger clustering than late- 
type (spiral) ones. Galaxy luminosity and color are 
perhaps the two major readily observed properties, 
which also facilitate comparison between different sur- 
veys and are less dependent on the stellar evolution 
models. Moreover, they have proven to be the two 
properties most predictive of galaxy environment, such 
that any residual dependence on morphology or sur- 
face brightness is weak (Blanton et al. 2005). In this 
paper, we measure the dependence of galaxy cluster- 
ing on color and luminosity for galaxies in the SDSS- 
III Baryon Oscillation Spectroscopic Survey (BOSS; 
Eisenstein et al. 2011; Dawson ct al. 2012) and study the 
implications. 

Galaxy clustering measurements, once properly inter- 
preted, can provide key information about galaxy for- 
mation and evolution. In particular, the theoretical un- 
derstanding of galaxy clustering has been greatly ad- 
vanced with the development of dark matter halo mod- 
els (see Cooray & Sheth 2002 and references therein). In 
the cosmological constant + cold dark matter (ACDM) 
paradigm, galaxies form and evolve in dark matter ha- 
los. While the properties and clustering of dark matter 
halos are well understood with the help of analytic mod- 
els and numerical simulations (e.g.. Mo & White 1996; 
Springel et al. 2005), the properties of galaxies are hard 
to predict because of the complex baryonic processes and 
the lack of complete theory of galaxy formation. 

Galaxy clustering offers an opportunity to connect 



galaxies to dark matter halos, providing a new direc- 
tion in studying galaxy formation and evolution. The 
halo occupation distribution (HOD) framework (see e.g.. 
Peacock & Smith 2000; Scljak 2000; Scoccimarro et al. 
2001; Berlind & Weinberg 2002; Berlind et al. 2003; 
Zheng et al. 2005) or the closely related conditional lumi- 
nosity function (CLF) method (Yang et al. 2003, 2005) 
describe the number of galaxies as a function of halo 
mass, and galaxy clustering is used to constrain the HOD 
or CLF parameters. The subhalo abundance matching 
(SHAM) method makes use of subhalos in high resolution 
iV-body simulations and connects them to galaxies to in- 
terpret galaxy clustering (see, e.g., Kravtsov et al. 2004; 
Conroy et al. 2006; Guo et al. 2010; Nuza et al. 2012). 
Such models essentially convert galaxy clustering mea- 
surements to the relation between galaxies and halos, 
which provides strong tests of galaxy formation models. 

The galaxy-halo connections inferred at different red- 
shifts, together with the theoretically known halo evo- 
lution, can lead to empirical constraints on galaxy evo- 
lution. For example, Zheng et al. (2007) compare the 
HODs for 0~1 DEEP2 and 2-0 SDSS galaxies and find 
a halo mass dependent growth of stellar mass of central 
galaxies, separating into contributions from star forma- 
tion and mergers. Tinker & Wetzel (2010) analyze four 
samples of galaxies from redshift z = 0.4 to z = 2.0. 
They find that more than 75% of the red satellite galax- 
ies move onto the red sequence because of halo mergers, 
while the mechanism for central galaxies to move to the 
red sequence evolves from z = 0.5 to z = 0. 

To advance our understanding of galaxy evolution, im- 
proved measurements of galaxy clustering in large galaxy 
surveys at different redshifts are necessary. The color and 
luminosity dependence of clustering has been studied in 
detail at z-^\ for DEEP2 galaxies (Coil et al. 2008) and 
at z^Q for SDSS galaxies (Zehavi et al. 2011; denoted 
Zll hereafter). The trends with color and luminosity are 
generally similar. However, while Coil et al. (2008) find 
no changes of the clustering within the red sequence, Zll 
find a continuous trend (in both amplitude and slope) of 
stronger clustering with color in SDSS galaxies. This 
discrepancy may be caused by different sample selec- 
tions, but it may also be a signature of galaxy evolution 
(e.g., related to the buildup of red sequence). Studying a 
sample from an intermediate redshift range, can provide 
new important information and a better understanding 
of galaxy evolution. 

In this paper, we use the recently released CMASS 
sample of the BOSS Data Release 9 (DR9; Ahn et al. 
2012) to measure the clustering of galaxies in the red- 
shift range 0.43 < z < 0.7 and study the dependence 
on galaxy luminosity and color. The sample is con- 
structed to contain a roughly volume-limited set of mas- 
sive and luminous galaxies (with a typical stellar mass 
of lO"-^/i"^M0; Marastonet al. 2012) in this redshift 
range. This sample has recently been used to accu- 
rately measure the baryon acoustic oscillation (BAO) 
signature (see Weinberg et al. 2012 for a recent com- 
prehensive review) on large scales (— 100 /i~^Mpc) by 
Anderson ct al. (2012). The sample has been thoroughly 
vetted, and the robustness of the results and cosmologi- 
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cal constraints from the BAO and redshift-space distor- 
tions are explored in a series of papers (Raid et al. 2012; 
Ross et al. 2012a,b; Sanchez et al. 2012; Samushia et al. 
2012; Scoccola et al. 2012; Tojeiro et al. 2012a,b). The 
smaller scale clustering measurements and HOD fits are 
first presented by White et al. (2011) (for an earlier 
smaller sample) and Nuza et al. (2012). Here we study 
the small to intermediate scale (0.05-40 /i~^Mpc) two- 
point correlation functions (2PCFs) of CMASS galax- 
ies, focusing especially on the dependence on luminosity 
and color and the implications on galaxy evolution from 
simple models. We will study the evolution of CMASS 
galaxies based on HOD modeling in a companion paper. 

The structure of the paper is organized as follows. In 
Section 2, we briefly describe the CMASS sample and 
our method of measuring the 2PCFs. The division to 
specific subsamples, the clustering measurements and de- 
tailed dependence on luminosity and color, and the im- 
plications for galaxy evolution are presented in Section 3. 
We summarize our results in Section 4. Appendix A dis- 
cusses the effect of different stellar evolution models on 
our results, and Appendix B explores the robustness of 
the jackknife error estimates used. 

Throughout the paper, we assume a spatially flat 
ACDM cosmology as in Anderson et al. (2012), with 
n,n = 0.274, h = 0.7, nth"^ = 0.0224, 7i, = 0.95, and 
tTg = 0-8, similar to the best-fit model from the WMAP 
7-year data (Komatsu et al. 2011). 

2. OBSERVATIONS AND METHODS 
2.1. Data 

As part of the SDSS-III survey, BOSS selects lu- 
minous galaxies from the multiple-band SDSS imaging 
(Fukugita et al. 1996; Gunn et al. 1998, 2006; York et al. 
2000) for spectroscopic observation to probe the large- 
scale BAO signals (Anderson et al. 2012). Dawson et al. 
(2012) provide a comprehensive overview of BOSS, while 
the technical details of BOSS are presented in Smee et al. 

(2012) and Bolton et al. (2012). The selection of BOSS 
galaxies is a union of targets in two different redshift in- 
tervals. One is an extension of the SDSS-I/II Luminous 
Red Galaxy (LRG) sample (Eisenstein et al. 2001), re- 
ferred to as LOWZ with 0.2 < z < 0.43 (Parejko et al. 
2012). The other, denoted as CMASS (White et al. 
2011; Anderson et al. 2012), includes --260,000 galax- 
ies in DR9, and is approximately stellar-mass limited at 
higher redshifts (0.43 < z < 0.7) with an effective volume 
of ~2.2Gpc'^ and an effective area of about 3, 300 deg^. In 
this paper, our study focuses on the CMASS sample. The 
target selection cuts for CMASS galaxies are fainter and 
bluer than the LRG sample in order to achieve a higher 
number density of about 3 x 10~^/i'^Mpc~'^ and sample a 
wider range of galaxies. The detailed target selection of 
the CMASS sample is described in Padmanabhan et al. 

(2013) , and summarized in Eisenstein et al. (2011) and 
Anderson et al. (2012). We briefly describe here the ma- 
jor selection cuts that will affect our analysis of luminos- 
ity and color dependence of the 2PCF. 

The CMASS sample aims at selecting galaxies fol- 
lowing a roughly constant stellar mass cut at redshift 
z > 0.4. The selection criteria of the CMASS galaxies 
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where d± is defined as 





d± = r„iod - «mod - (5mod - fmod)/8 (6) 

and all magnitudes are extinction corrected 
(Schlegel et al. 1998) and are in the observed frame. 
While the magnitudes are calculated using cmodel 
magnitudes (denoted by the subscript "cmod"), the 
colors are computed using model magnitudes (denoted 
by the subscript "mod" ; see Anderson et al. 2012 for 
details). The magnitude ifib2 corresponds to the z-band 
flux within the 2" fiber size. CMASS objects also pass 
specific star-galaxy separation cuts, as described in 
Anderson et al. (2012). 

2.2. Selection Cuts and Sample Completeness 

The CMASS galaxies are chosen by applying complex 
target selection cuts, and it is thus difficult to construct 
exact volume-limited samples. Since we intend to in- 
vestigate the luminosity and color dependence of galaxy 
clustering in this paper, we must pay particular atten- 
tion to the completeness in luminosity and color. With 
this goal in mind, we first investigate the target selection 
cuts projected to the color-luminosity plane at different 
redshifts, which will provide the appropriate boundaries 
in constructing our samples. 

Figure 1 shows the color-magnitude diagram (CMD) 
at six narrow redshift ranges with Az ~ 0.05. The 
contours represent the number density distribution of 
CMASS galaxies in the CMD. The absolute magnitude 
Mi and r — i color are k-|-e corrected to z = 0.55 
throughout the paper using a global Flexible Stellar 
Population Synthesis (FSPS) model (Conroy et al. 2009; 
Conroy & Gunn 2010; Tojeiro et al. 2011, 2012b). Using 
other stellar evolution models slightly changes the result- 
ing magnitudes and colors, but does not affect our anal- 
ysis of galaxy clustering, as we discuss in Appendix A. 
In each panel of Figure 1, the approximate positions of 
the selection cuts are marked as dashed lines. We fo- 
cus on the three main cuts (Equations 1-3). Following 
Zehavi et al. (2011), we adopt a luminosity-dependent 
color cut to separate red and blue galaxies (discussed in 
detail in Section 3.3), 

(r - i)cut = 0.679 - 0.082(A/i + 20), (7) 

represented by the green lines in Figure 1. 

The horizontal cut in each panel corresponds to the 
z-band faint-end flux limit (Equation 1, z < 19.9), which 
selects galaxies brighter than the corresponding abso- 
lute magnitude at each redshift. The slightly tilted 
vertical dashed line on the left of the distribution is 
the d± cut (Equation 2), which removes galaxies with 
bluer colors and with lower redshifts (Cannon et al. 2006; 
Padmanabhan et al. 2013). The bottom-left dashed line 
is the z-band sliding cut (Equation 3), which excludes the 
fainter and bluer (thus lower stellar mass) galaxies from 
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Fig. 1. — Color-magnitude diagram (CMD) of CMASS galaxies in different redshift intervals. Both magnitude and color arc k+c corrected 
to z = 0.55. The contours represent the number density distribution of CMASS galaxies in the CMD. The approximate positions of the 
selection cuts are shown as dashed lines in each panel, labeled with the corresponding equation numbers of the selection cuts (see the text). 
The green solid lines are our color cut (Equation 7) for red and blue galaxies. (A color version of this figure is available in the online 
journal.) 



the sample. These three main cuts evolve with redshift. 
A consequence of the d± and sliding cuts is that there 
are more blue galaxies at higher redshifts. At z > 0.55, 
while the CMASS sample is less affected by the d± cut, 
the sliding cut must be carefully taken into account at 
all redshifts in constructing complete galaxy samples. At 
z < 0.55, blue galaxies are highly incomplete (see an esti- 
mation of the fraction of star-forming galaxies in CMASS 
in Figure 16 of Chen et al. 2012). At z < 0.45, the di_ 
cut leads to incompleteness even for the most luminous 
red galaxies. 

A more sophisticated method to study the sample 
completeness would be to simulate the galaxy proper- 
ties (stellar mass, luminosity and color) at different red- 
shifts by assuming certain galaxy stellar evolution models 
and to apply the the selection cuts to simulated galaxies 
(Swanson et al. 2013). While it would generally depend 
on the assumptions in the galaxy evolution models, such 
a study can provide a quantitative estimate of the sam- 
ple completeness as a function of color and luminosity. 
In this paper, however, we do not intend to evoke such a 
model in our clustering analysis. Instead, we empirically 
use the CMD and selection cuts at each redshift to con- 
struct approximately complete galaxy samples. As we 
proceed in our analysis, we keep in mind the boundaries 
of "completeness" defined by the selection cuts when con- 
structing our galaxy samples. One advantage of such an 
empirical method is that it is largely model- independent. 

2.3. Clustering Measurements 



In this paper, we focus our discussion on the 
galaxy 2PCF. We use the Landy-Szalay estimator 
(Landy & Szalay 1993) to measure the 2PCF of galax- 
ies, 

^, ^ DD-2DR + RR 

= RR 

where DD, DR and RR are the data-data, data-random, 
and random-random pair counts measured from the data 
of N galaxies and random samples consisting of Nn ran- 
dom points. These pair counts are appropriately normal- 
ized by N{N - l)/2, NNr, and Nii{Nr - l)/2, respec- 
tively. 

We measure the 3-dimensional (3D) 2PCF ^(r-pjT^), 
where rp and r^r are the separations of galaxy pairs per- 
pendicular and parallel to the line of sight. The redshift- 
space 2PCF ^(r-p,r^) differs from the real-space one be- 
cause of redshift distortions induced by galaxy peculiar 
velocities. The redshift distortions can be mitigated by 
projecting the 2PCF along the line-of-sight direction, 
with the projected 2PCF Wp{rp) (Davis & Peebles 1983) 
defined and measured as 

/•oo 

Wp{rp)^2 / ^(rp,r^)dr^ = 2 V'$(rp,r^.i)Ar^_j(9) 
Jo 

where r7r,i and Ar-^^ are the i-th bin of the line-of-sight 
separation and its corresponding bin size. In practice, 
wc sum ^{rp,r-^) along the line-of-sight direction up to 
?'7r,max = 80 /i~^Mpc to iucludc most of the correlated 
pairs. As we measure the 2PCFs up to rp ~ 40 /i~^Mpc, 
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the clustering measurements do not depend significantly 
on the assumed rTr^max once it is sufficiently larger than 
40 /i~^Mpc. For example, if we integrate the line-of- 
sight direction to 200 /i^^Mpc, the contribution from 
80/i"^Mpc < < 200/i-^Mpc to Wp{rp) only intro- 
duces noisy fluctuations of about 2%. 

The projected 2PCF can be related to the real-space 
correlation function, ^(r), by 



Wp{rp) 



rdr£,{r){r'^ - r^) 



-1/2 



(10) 



(Davis & Peebles 1983). It is common practice to char- 
acterize ^(r) by a power law on small scales, ^(r) = 
i'r/ro)~'^, and in such a case Wp{rp) can be expressed 
as 



Wp{rp) 



^l'rri)r(l, 



Til) (11) 



The covariance matrix of the correlation function is es- 
timated using the jackknife resampling method (follow- 
ing Zehavi et al. 2005b, 2011) 



cov(6,e,) = ^T.(^l~mh^~^) 



N 



N 



(12) 



1=1 



where N = 100 is the total number of jackknife samples, 
is the 2PCF in the i-th pair separation bin measured 
from the l-th jackknife sample, and £,i is the average 
over all samples. We expand on the tests performed by 
Zehavi et al. (2005b) and investigate the accuracy of the 
jackknife error estimates using CMASS mock catalogs 
and find good agreement between the jackknife estimates 
and those from multiple mocks (see Appendix B for more 
details) . 

We focus our measurement of the galaxy 2PCF on 
the small to intermediate scales (0.05 h~-^Mpc < rp < 
40 /i~"'^Mpc), aiming at probing the relation between 
galaxies and their host halos. However, the small-scale 
measurement of the 2PCF is limited by the fiber col- 
lision effect in the spectrograph system of SDSS-III, 
where two fibers on the same plate cannot be placed 
closer than an angular separation of 62" (Dawson et al. 
2012; Anderson et al. 2012). As a result, about 5.5% 
of CMASS galaxies do not have redshifts, strongly af- 
fecting small-scale clustering measurements. We correct 
for this effect using the method proposed and tested by 
Guo et al. (2012), which divides the galaxy sample into 
two distinct populations, one free of fiber collisions (re- 
ferred to as Di) and the other consisting of potentially 
collided galaxies (denoted as D2). The total clustering 
signal is a combination of the contributions from these 
two populations, where the contribution of the collided 
population is estimated from the resolved galaxies in tile- 
overlap regions. 

As discussed in detail in Guo et al. (2012), there are 
two main systcmatics that could impact the accuracy of 
such a method to treat fiber collisions. One is possi- 
ble density variations between the tile overlap and non- 
overlap regions, which is found to be insignificant from 
tests with mock galaxy catalogs. Another effect is that 
galaxies in collided triplets (or even higher-order collid- 
ing groups) can only be fully recovered in regions cov- 



ered by at least three tiles, making the estimation from 
close pairs in two-tile regions not accurate (because of 
the lack of D2D2 close pairs). This effect is alleviated by 
an additional correction term using the measured close 
pairs in the two-tile regions. After full application of the 
Guo et al. (2012) method, we estimate the remaining sys- 
tematic errors to be less than 3-5%. It is generally diffi- 
cult to have an unbiased correction on all scales, but this 
approach provides the best estimate of the galaxy 2PCF 
on small scales, compared with other possible methods. 

When counting pairs for the 2PCF, each galaxy is 
assigned a series of weights to reduce variance in the 
measurements and take into account different effects. 
Following Anderson et al. (2012), we apply a scale- 
independent weight to optimize the clustering measure- 
ments (Feldman et al. 1994) 



1 



wfkp 



1 + n{z)Po ' 



(13) 



where n{z) is the mean density at redshift z, and Pq = 
2 X 10^/i~'^Mpc'^. This equation provides similar results 
to the minimum variance weight used by Hamilton (1993) 
and Zehavi et al. (2002). Another employed weight, Wif , 
accounts for the fact that not every galaxy with a spec- 
trum taken has a reliable redshift measurement. The 
"redshift failures" are dependent on the positions of 
fibers on the plates (Ross et al. 2012b) and are corrected 
by up-weighting the nearest galaxy that has an accu- 
rate redshift. The final weight, Wgys, is caused by the 
scarcity of galaxies detected due to foreground bright 
stars (Ross et al. 2011b). Ross et al. (2012b) present 
a comprehensive study of potential systematic effects 
in the 2PCF analysis, and compute a set of weights 
Wsys based on the stellar density and ifib2 magnitude. 
Therefore, the total weight applied to each galaxy is 



Wtot = WpKPWsysWrf- 



(14) 



The quantity wtot is applied to both the Di and D2 pop- 
ulations in the fiber collision correction. The systematic 
weight Wsys only has a small effect on small and inter- 
mediate scales, but significantly changes the clustering 
measurements on BAG scales. 

We construct the random catalogs according to the de- 
tailed angular selection of the DR9 galaxy sample. The 
radial selection function for each sample is taken into ac- 
count by assigning the shuffled galaxy redshifts to the 
random objects. The shuffling method provides a better 
representation of the true distribution compared with a 
smooth spline fit to the observed galaxy redshift distri- 
bution, as detailed in Ross et al. (2012b). This process is 
done separately for the northern and southern Galactic 
Caps to account for the different number density distri- 
butions (see Anderson et al. 2012, for details). To apply 
the fiber collision correction, separate random catalogs 
for the Di and D2 populations are used. Denoting the 
fraction of recovered D2 galaxies as N2/N2, for the D2 
random catalogs we apply an additional angular mask 
N2/N2 in each sector (see Guo et al. 2012 for details). 

3. RESULTS 

3.1. 2PCF of the Full CMASS sample 

Before presenting the luminosity and color dependence 
of the 2PCFs for CMASS galaxies, we show in Figure 2 
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Fig. 2. — Projected (left panel) and redshift-spacc (right panel) 2PCFs for the entire CMASS sample in the range of 0.43 < z < 0.7. The 
solid lines present our measurements. The open symbols are the recent measurements of Nuza et al. (2012), which are in good agreement. 
The vertical dashed lines correspond to the maximal fiber collision scale (~0.53 /i~^Mpc) for z = 0.7. The dotted line in the left panel is a 
power-law fit to Wp for 0.05/i~^Mpc < rp < 25/i~^Mpc with the corresponding parameters as labeled. 



the projected and redshift-space 2PCFs for the entire 
CMASS sample in the redshift range 0.43 < z < 0.7 
(soHd hnes). For comparison, we also display the re- 
cent CMASS DR9 measurements of Nuza et al. (2012, 
open symbols), which are limited to slightly larger scales 
(> 0.5/i~^Mpc) and show good agreement on all mea- 
sured scales. The vertical dashed lines indicate the max- 
imal fiber collision scale, ^^0.53 /i~^Mpc at redshift z = 
0.7. By applying the fiber-collision correction method of 
Guo et al. (2012), we are able to robustly measure the 
small-scale clustering (note the small error bars). This 
should enable better constraints on the spatial distribu- 
tion of galaxies inside dark matter halos and a better 
constraint on the fraction of satellite galaxies, which we 
will address in future work. 

The dotted line in the left panel of Figure 2 is the 
power-law fit to Wp for 0.05 h~^Mpc < r^, < 25 ft.~^Mpc, 
based on Equation 11, using the full error covariance 
matrix. Although Wp clearly deviates from a power-law 
(also shown from the x^/dof of the fitting), a power- 
law fit provides a simple characterization of the cluster- 
ing and allows for easy comparisons among the 2PCFs 
of different galaxy samples. The correlation length for 
the CMASS sample is ro = 7.93 ± OMh-^Mpc and the 
slope 7 = 1.85 ± 0.01 (note that given the large best- 
fit x^, the error bars here should be interpreted with 
care). The amplitude is slightly larger than that of the 
Mr < —21.5 luminosity-threshold SDSS DR7 main sam- 
ple (Zll). This is consistent with the expectation that 
the CMASS galaxies occupy the bright end of the lumi- 
nosity function. These values arc also similar to cluster- 
ing measurements of LRGs at z ^ 0.55 from the 2dF- 
SDSS LRG and QSO (2SLAQ) survey (Cannon et al. 
2006; Ross et al. 2007; Wake et al. 2008), as expected, 
since the CMASS galaxy color selections were based on 
the 2SLAQ LRG selection. 

By fitting the ratio between the measured galaxy Wp 
and the theoretically computed dark matter Wp on scales 



Sh-^Mpc < rp < 25/i"iMpc (detailed in Section 3.5), 
we obtain a "large-scale" galaxy bias factor of 6 = 2.16± 
0.01, consistent with the measurements of other authors 
(White et al. 2011; Nuza et al. 2012; Sanchez et al. 2012; 
Shen et al. 2012). The galaxy bias determined from the 
projected 2PCF has less scale dependence than that from 
the redshift-space 2PCF. The exact value of the implied 
bias can depend on the fitting scales and fitting meth- 
ods. We find that if the minimal fitting scale is changed 
from 3/i~^Mpc to 5/i~^Mpc (or larger), the resulting 
bias only varies slightly at the 2a level. In the right 
panel of Figure 2 we present the redshift-space 2PCF 
^(s) on scales above s > 0.1 /i~^Mpc, where it can be re- 
liably measured. The fiber collisions in that case impact 
larger scales than indicated by the dashed line, since s 
includes the contribution from the line-of-sight separa- 
tions, = -t- (see more discussion in Guo et al. 
2012). 

3.2. Luminosity Dependence 
3.2.1. Luminosity Cuts 

We now investigate the luminosity dependence of 
CMASS galaxy clustering. To minimize the influence of 
sample incompleteness, we carefully construct samples 
of different luminosities by accounting for the selection 
cuts as a function of redshifts discussed in Section 2.2. 
Wc divide galaxies into two redshift bins, 0.43 < z < 0.55 
and 0.55 < z < 0.7. The color- magnitude distributions 
in these two redshift bins, the overall magnitude-redshift 
distribution, and the cuts used to define our samples are 
shown in Figure 3. In the right-most panel, some galaxies 
lie below the faint fiux limit (denoted by the lower yellow 
line), reflecting the change between the photometry at 
the time of targeting and that from the final processing. 
To keep a uniform criterion, we construct our luminosity 
samples based on the targeting photometry (see details 
in Anderson et al. 2012). We avoid the impact of the 
sliding cut by only selecting galaxies brighter than the 
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Fig. 3. — Color-magnitude diagrams of CMASS galaxies in the two redshift bins we use, as well as overall distribution of galaxies in 
i-band absolute magnitude and redshift. The red lines delineate the luminosity bin samples we study. The two green lines in the right 
panel represent the j-band flux limits of Equation 1, which is also k+e corrected to z = 0.55. (A color version of this flgure is available in 
the online journal.) 



intersections between the d± and sliding cuts. The in- 
completeness caused by the d± cut at z < 0.55 cannot be 
avoided. Such a limitation means that the blue galaxies 
are incomplete, especially for the low-redshift samples, 
while the red galaxies are close to complete for z > 0.5, 
which is a caveat to remember when interpreting the re- 
sults. We will thus also study the luminosity dependence 
limited to the more complete red galaxies. The sliding 
cut also impacts our ability to study fainter galaxies at 
high redshift, resulting in the unsampled "triangle" re- 
gion above z = 0.55 in the right panel of Figure 3. We 
have three and two luminosity bin samples at lower and 
higher redshifts, respectively, each with a bin width of 
0.3 magnitude, as shown in the figure. 

The total numbers of galaxies, iVtot, and the comov- 
ing volume, Vz, in each luminosity bin are shown in 
Table 1. We also provide the number of blue and 
red galaxies using the color cut (Equation 7). Power- 
law fits of Equation 11 to the projected 2PCF Wp for 
0.1/i~^Mpc < Tp < 25/i~^Mpc arc also given in the ta- 
ble. It is evident that at lower redshifts, red galaxies 
dominate the samples, and blue galaxies contribute less 
than 10%. At higher redshifts, approximately 32% of 
the sample are blue galaxies. This could still be partly 
related to the CMASS selection effects. If the change 
in blue galaxy fraction is caused by evolution, i.e., blue 
galaxies turning red with time, we would expect the lu- 
minosity dependence of the 2PCF for red galaxies to also 
evolve with redshift. 

3.2.2. The Dependence of Galaxy 2PCF on Luminosity 

The projected 2PCFs of different luminosity samples 
at the two redshift bins are shown in the top panels of 
Figure 4. At both redshifts, the luminosity dependence 
of Wp{rp) is evident, with more luminous galaxies ex- 
hibiting stronger clustering, consistent with the results 
of the SDSS-I/II main sample (Zehavi et al. 2005b, 2011; 
Li et al. 2006). At 2 < 0.55, since red galaxies contribute 
90% of the CMASS galaxy population, the luminosity 
dependence mostly reflects the clustering environment of 
the red galaxies. At z > 0.55, our measurements of the 
2PCF become noisier because of the lower number of 
galaxies. Even after accounting for the uncertainties in 



the measurements, the luminosity dependence of cluster- 
ing is significant in the redshift range of CMASS galaxies. 

According to the HOD modeling results in SDSS-I/II 
(Zehavi et al. 2005b, 2011), the increase in the clustering 
amplitude for more luminous samples reflects the shift in 
the host halo mass towards the high mass end. The mod- 
eling results in Zehavi et al. (2005b, 2011) also show that 
the satellite fraction drops as the luminosity of galaxies 
increases. Our measurements naively seem to support 
such a result - although the 2PCFs becomes generally 
noisier for higher luminosity samples, the uncertainties 
in the 2PCFs on small scales increase faster. At such 
scales, satellite galaxies have a significant contribution 
to the clustering signal. Therefore, the increase in the 
measurement errors could be a reflection of the lack of 
satellites in more luminous samples. 

The shapes of Wp for all luminosity samples are simi- 
lar. The deviation from a power-law in Wp is somewhat 
more apparent for brighter galaxies, consistent with the 
results from main galaxies (Zehavi et al. 2004, 2005b, 
2011). In the halo model, the slope of Wp has a rapid 
change around a few Mpc, indicating the transition from 
intra-halo galaxy pairs (1-halo term) to inter- halo galaxy 
pairs (2-halo term). In our measurements, we see that 
this transition scale increases with increasing luminos- 
ity, in agreement with the interpretation that more lu- 
minous galaxies reside in more massive (hence larger) 
halos. There is a weak trend that brighter galaxies have 
a steeper slope in Wp on small scales, in line with the re- 
sult in Zll. The trend can also be interpreted as a result 
of the change in the host halo mass scale (see Fig. 7 and 
Appendix A in Zheng et al. 2009). 

The bottom panels of Figure 4 show the redshift evo- 
lution of the 2PCFs of galaxies in a fixed luminosity bin 
-22.5 < Afj < -22.2. The left panel is for all the galax- 
ies (both blue and red). Galaxies at a lower redshift 
appear to have a higher clustering strength. This re- 
sult may be due to the incompleteness of blue galaxies 
at lower redshifts, and the inclusion of more (less clus- 
tered) blue galaxies at higher redshifts. We therefore 
also examine the redshift evolution of the 2PCF in the 
—22.5 < Mi < —22.2 sample for red galaxies only, as 
shown in the bottom-right panel of Figure 4, where the 
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TABLE 1 

Samples of different luminosities 



Mi range 


z range 








y,[Gpc/h]^ 


ro 


7 


X^/dof 


-21.9, 


-21.6 


0.43, 0.55 


48391 


46235 


2156 


0.484 


7.64 ±0.12 


1.86 ±0.02 


47.89/10 


-22.2, 


-21.9 


0.43, 0.55 


24190 


22419 


1771 


0.484 


8.49 ±0.18 


1.91 ±0.03 


45.87/10 


-22.5, 


-22.2 


0.43, 0.55 


9308 


8687 


621 


0.484 


9.99 ±0.31 


1.89 ±0.04 


10.22/10 


-22.5, 


-22.2 


0.55, 0.70 


23404 


15821 


7583 


0.851 


8.56 ±0.19 


1.91 ±0.03 


10.68/10 


-22.8, 


-22.5 


0.55, 0.70 


7484 


5135 


2349 


0.851 


10.40 ±0.32 


1.88 ±0.06 


15.13/10 



Note. — ro and 7 are obtained from fitting a power-law to Wp{rp) using the full error covariance 
matrices for 0.1 h~^lS/lpc < rp < 25 /i~^Mpc. The ratios between 3,nd degrees-of-freedom (dof) of the 
fits are also shown. 
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Fig. 4. — Projected correlation functions, Wp{rp), for the various luminosity subsamples at low (top left) and high rcdshift (top right). 
The bottom panels present the redshift evolution oiwp{rp) in the luminosity interval —22.8 < Mi < —22.2 for all the galaxies in the sample 
(left) and only for the red galaxies (right). Error bars shown are from the diagonal elements of the jackknife covariance matrices. (A color 
version of this figure is available in the online journal.) 
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red galaxies are defined by the color cut in Equation 7. 
We find only slight evolution with redshift for the red 
galaxies (at most 18% in Wp for Tp > l/i~^Mpc), not 
significant within the measurement errors, implying that 
the differences in the sample of all galaxies (bottom-left 
panel) are mostly induced by the blue galaxies. 

In the redshift range of 0.16 < z < 0.44, Zehavi et al. 
(2005a) find no strong evolution trend in Wp{rp) in the 
SDSS LRG sample. Their sample of -23.2 < Mg < 
—21.8 at 0.16 < z < 0.44 has a comoving number density 
of n~0.2 X lO-^'/i^Mpc"^ (their Figure 2), similar to the 
number density of our sample of —22.5 < Mi < —22.2 
at 0.43 < z < 0.7 shown in the bottom right panel of 
Figure 4. Combining their results with ours, we would 
infer that there is no strong redshift evolution in Wp of lu- 
minous red galaxies in the redshift range of 0.1 < z < 0.7, 
consistent with the results of Wake et al. (2008), imply- 
ing that the effect of structure growth roughly cancels 
that of evolution of galaxy bias. As will be discussed 
in Section 3.5, within the errorbars, the trend is also 
roughly consistent with passive evolution. 

3.3. Color Dependence 
3.3.1. Color Cuts 

The various target selection cuts make it difficult to 
discern the "red sequence" and "blue cloud" in the color- 
magnitude diagram of Figure 1. Motivated by the bi- 
modal color distribution of galaxies (e.g., Strateva et al. 
2001; Baldry et al. 2004), we construct the division be- 
tween red and blue galaxies by fitting the galaxy color 
distribution with two Gaussian distributions. Figure 5 
shows the probability distribution function of r — z color 
(k+c corrected to z=0.55) in six redshift slices as in 
Figure 1. The distribution at each redshift is computed 
from galaxy samples that are as complete in luminosity 
as possible, i.e., we only use galaxies more luminous than 
the luminosity given by the intersection of line (2) and 
line (3) in Figure 1, which corresponds to an z-band ap- 
parent magnitude of icmod = 19.46. The CMASS sample 
shows a clear bimodal distribution in color, similar to the 
findings in SDSS-I/II (Strateva et al. 2001; Baldry et al. 
2004; Skibba & Sheth 2009). We can use the intersec- 
tion of the best-fit two Gaussian distributions to divide 
galaxies into blue and red samples. 

We do not perform the two-Gaussian fit to the galaxy 
distribution in redshift interval 0.4 < z < 0.45, as these 
galaxy samples are far from complete. The blue galax- 
ies are essentially missing from this sample, and even 
the red galaxy colors at this redshift are not well de- 
scribed by a Gaussian distribution. As shown in Figure 1. 
galaxies at this redshift suffer from the d± selection cut, 
which eliminates the blue galaxies and a fraction of red 
ones. In the 0.45 < z < 0.5 redshift bin, we still miss 
galaxies with r — i < 0.9, and the distribution is domi- 
nated by the Gaussian profile from the red population. 
At 0.5 < z < 0.55, luminous blue galaxies are excluded 
from the sample by the d± cut (see Figure 1), and the 
contribution to the blue Gaussian profile is mainly from 
faint blue galaxies. Therefore, at z < 0.55, blue galax- 
ies in the GMASS sample are far from complete. For red 
galaxies, we find that the centers of the color distribution 
do not significantly change with redshift. Thus for the 
analysis of the whole CMASS sample, we use the redshift- 



independent color cut for simplicity (see Equation 7). 

The red/blue color division cut shows a mild redshift 
dependence, becoming bluer at higher redshift. Since 
both the color and magnitude used in this paper have 
been k-fe corrected (i.e., the evolution effects are re- 
moved), such a mild evolution might indicate that the 
global evolution correction is not accurate. On the other 
hand, the photometric errors increase for larger redshift 
(see below), making the two Gaussian profiles broader, 
which can lead to a shift of the red-blue division cut 
towards the blue end even if there is no change in the 
blue and red populations. Moreover, the blue sample is 
generally incomplete due to the selection effects, which 
may also introduce additional change of the color cut. 
Therefore, the weak dependence of the red-blue division 
cut on redshift may not reveal much about the evolution. 

We further investigate the dependence of the color dis- 
tribution on luminosity for the whole CMASS sample 
at 0.43 < z < 0.7. As shown in Figure 6, the peak 
of the red sequence, as well as the intersection of the 
two Gaussian profiles, become slightly redder as the lu- 
minosity increases, refiecting the well-known tilt of the 
red sequence. The tilted red sequence is likely a refiec- 
tion of differences in the chemical composition, where 
the more luminous galaxies are richer in metals while the 
smaller galaxies suffer from the loss of metal-enriched gas 
(Kodama & Arimoto 1997; Gallazzi et al. 2006). The 
tilted red sequence may also reflect the role of dry merg- 
ers (i.e., of gas-poor galaxies) in the evolution of red 
galaxies (e.g., Skelton et al. 2009), with galaxies increas- 
ing their mass (luminosity) from mergers and becom- 
ing older (redder) as a result of stellar evolution (see 
Faber et al. 2007 for a comprehensive review). 

In order to find a reasonable color cut for red and blue 
galaxies in CMASS, we fit the bimodal color distribution 
as a function of luminosity for galaxies in the range of 
0.5 < z < 0.7, where the samples are less affected by in- 
completeness. The resulting color cut is the one already 
presented in Equation 7. With such a color cut, if we 
naively count the CMASS galaxies disregarding the in- 
completeness at lower redshifts, we find that only about 
13% of the galaxies in CMASS are blue galaxies, and 80% 
of these blue galaxies are from z > 0.55. We emphasize 
that our color cut is based on the k-l-e corrected colors. 

Masters et al. (2011) proposed an observer-frame color 
cut oi g — i = 2.35 for CMASS galaxies, motivated by the 
color and morphology distribution for a matched sample 
between CMASS and Hubble Space Telescope imaging of 
the Cosmic Evolution Survey (COSMOS; Scoville et al. 
2007). They demonstrated that the g — i cut can be used 
to separate the elliptical and spiral galaxies in CMASS. 
However, there is still some fraction of LRG progenitors 
with g — i < 2.35 due to small amounts of star forma- 
tion (Tojeiro et al. 2012b). For our purpose of study- 
ing the color dependence of clustering, an observer-frame 
color cut would mix galaxies from different populations, 
since galaxies of the same intrinsic colors will appear to 
have different observed colors at different redshifts. In 
addition, the g-band magnitude in CMASS is especially 
faint and thus large measurement errors will make g-band 
based colors prone to spurious fluctuations. We therefore 
prefer to use the "intrinsic" (i.e., k+e corrected) r ~ i 
color in our study. We note that the z = 0.55 observed 
r- and i-bands are close to rest-frame g- and r-bands, 
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Fig. 5. — Probability distribution function of r — i color at different redshift slices, for CMASS samples of high completeness (see text). 
The black lines are the histogram of r — i for all galaxies. The red and blue lines are the bimodality fitting using two Gaussian distributions, 
with the green curves as their combination. We do not fit the distribution of r — i for 0.4 < z < 0.45 because both the red and blue galaxies 
are far from complete in this redshift interval. (A color version of this figure is available in the online journal.) 
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Fig. 6. — Similar to Figure 5, but for different luminosity intervals at 0.43 < z < 0.7. (A color version of this figure is available in the 
online journal.) 

respectively. So our analysis based on z = 0.55 Mi and 
r — i approximates that of rest-frame Mr and g — r, as 
adopted in Zll for 2; ~ SDSS galaxies. 

To study the color dependence in different luminos- 
ity and redshift intervals in more detail, we further de- 
compose the sample into finer color subsamples. The 
Gaussian fittings provide the centers and la widths of 
the blue cloud and red sequence, which are used in defin- 
ing the fine color cuts 



With the three cuts, we can form blue (below the be cut), 
green (between the be and br cuts), redseq (between the 
br and rr cuts), and reddest (above the rr cut) samples. 
In each redshift interval, the luminosity-dependent color 
cuts are fitted with a straight line. 



r — i = ajMi + bj 



(18) 



(r — i)bc = blue center (15) 
{r — i)hi- = I'cd center — 0.5 x red width (16) 



where j = 1, 2, 3 for the 6c, br, and rr cuts, respectively. 
The linear fits for these cuts are listed in Table 2. For 
clarity, we show again the color-magnitude diagram in 
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TABLE 2 

- i COLOR CUTS FOR DIFFERENT REDSHIFT INTERVALS 



2 range 


0,1 


bi 


0.2 


b2 




63 




^ green 


-^rcdscq 


-^reddest 


y4Gpc/h]3 


0.45, 0.50 


-0.023 


0.386 


-0.0199 


0.557 


-0.0089 


0.849 


934 


13233 


20047 


25094 


0.19 


0.50, 0.55 


-0.099 


-1.315 


-0.0622 


-0.398 


-0.0346 


0.281 


6215 


17415 


24721 


27361 


0.22 


0.55, 0.60 


-0.249 


-4.777 


-0.057 


-0.294 


-0.0251 


0.496 


2841 


20323 


19520 


18468 


0.25 


0.60, 0.65 


-0.107 


-1.619 


-0.054 


-0.227 


-0.0242 


0.538 


3892 


13349 


12261 


9374 


0.28 


0.65, 0.70 


-0.034 


-0.074 


-0.0745 


-0.750 


-0.0547 


-0.193 


2517 


6517 


5797 


5236 


0.31 



Figure 7, with the fine color cuts superimposed. The 
redshift interval of 0.4 < z < 0.45 is omitted because of 
the high sample incompleteness. 

Examining Figures 5-7, there is a noticeable trend that 
the width of the red sequence appears narrower for more 
luminous galaxies and at lower redshifts. We highlight 
this effect in Figure 8, which presents the Icr width of the 
red sequence as a function of luminosity and redshift. 
Since the photometric errors become larger for fainter 
galaxies and at higher redshifts, we subtract their contri- 
bution in quadrature to obtain the intrinsic color width 
of the red sequence galaxies. The r — i color photometric 
errors are estimated by simply combining in quadrature 
the errors in r- and i-band magnitudes (neglecting any 
additional errors in the k-l-e corrections and in the cor- 
relation between r- and z-band photometric errors). As 
shown in Figure 8, the above trend persists for the intrin- 
sic color scatter, suggesting an evolutionary effect which 
we discuss further in Appendix A. 

3.3.2. The Dependence of Galaxy 2PCF on Color 

With the color cuts defined in the previous subsec- 
tion, we investigate the dependence of galaxy 2PCFs on 
the r — i color. First, we examine the 2PCFs for blue 
and red galaxies in the whole CMASS sample. The red 
and blue samples here are defined using the color cut 
in Equation 7. The samples are fiux-limited (in addi- 
tion to the other selection cuts) and are by no means 
complete. The purpose of this exercise is simply to have 
an overall view of the difference in red and blue galaxy 
clustering. The ^(''pj^ir) measurements for blue and red 
galaxies are shown in Figure 9. Red galaxies are more 
strongly clustered. The "Finger-of-God" feature on small 
scales caused by random motions of galaxies in virial- 
ized structures can be clearly seen for both red and blue 
galaxies. Red galaxies have a stronger "finger-of-God" 
effect, reflecting their stronger motions within halos. On 
large scales (e.g., above rp = 10/i~-^Mpc, the outmost 
contours), the contours for both blue and red galaxies 
show the flattening trend caused by coherent large-scale 
infall (Kaiser 1987). On these scales, the Kaiser squash- 
ing effect appears to be stronger for blue galaxies, since 
the effect is determined by w fl'^^ /b and blue galaxies 
have a smaller galaxy bias factor b. 

We now investigate the color-dependent 2PCFs as a 
function of luminosity and redshift from the fine color 
samples. In order to minimize the effect of incomplete- 
ness, the luminosity and redshift bins are selected using 
Figure 1 to make sure that red galaxies are not affected 
by the selection cuts. The blue galaxies are generally not 
complete at most redshifts and the results of the blue 
samples need to be interpreted with care. Nevertheless, 
the blue samples are still useful in comparison with the 
red galaxies. 



The main results of the color-dependent 2PCFs are 
summarized in Figure 10. The top panels display the de- 
pendence of Wp on color in the magnitude range —22.5 < 
Ali < —21.5 at two different redshift intervals. The trend 
with color is obvious at both redshifts - there is a contin- 
uous increase in the clustering amplitude as galaxy color 
goes from blue to red. This result is consistent with the 
behavior observed in the SDSS-I/II main galaxy sam- 
ple (Zll). On small scales (below the inflection scale of 
l-2/i~^Mpc), there appears to be a trend that redder 
galaxies have a steeper slope in Wp, which is weaker than 
that measured by Zll. According to the HOD modeling 
result in Zll, for galaxies in a flxed luminosity range, red- 
der galaxies generally have a higher fraction of satellites 
residing in massive halos. Our results therefore implies 
that a larger fraction of redder galaxies reside in more 
massive halos, giving rise to a larger clustering ampli- 
tude. The steepening of Wp on smaller scales may also 
indicate a halo mass scale shift with color, leading to 
a relative increase in the contribution from the 1-lialo 
central-satellite pairs with respect to the 1-halo satellite- 
satellite pairs (see the Appendix A of Zheng et al. 2009). 

The bottom panels of Figure 10 present the 2PCFs 
for more luminous red galaxies, with —23 < Mi < —22, 
showing the results for the redseq and reddest color sub- 
sets. The number densities of these two subsets have 
a roughly constant value of 0.21 x 10~'^h^Mpc~'^ and 
0.18 X lO^^ft.^Mpc"'^ over the redshift range, respectively. 
There is a stronger inflection in the slope of Wp below 
2ft,~^Mpc for these luminous galaxies and a trend of 
steeper slope for the reddest sample than for the red- 
seq sample, indicating the shift in the halo mass scale. 
The two bottom panels also compare the redshift evolu- 
tion for galaxies of the same color and luminosity. No 
strong evolution in Wp is found, consistent with the bot- 
tom right panel of Figure 4. The implications of these 
results are discussed further in Section 4. 

3.4. Red Galaxy Samples with Fixed Number Density 

In previous sections, we constructed galaxy samples in 
certain luminosity and color bins, in an attempt to min- 
imize the influence of incompleteness caused by target 
selection cuts. Motivated by a simple passive evolution 
model, we can further define galaxy samples with fixed 
number density at different redshifts (White et al. 2007; 
Brown et al. 2008; Wake et al. 2008). If during the evo- 
lution each galaxy in a sample retains its identity, ex- 
periencing no merger or disruption, the number density 
of the galaxy sample would not change with redshift. 
The evolution of the 2PCF of such a galaxy sample can 
be readily predicted (Fry 1996). If, in addition, no star 
formation occurs in these galaxies during the process, 
their stellar population would evolve passively and can be 
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Fig. 7. — Our adopted color cuts of Table 2 in the color-magnitude diagram. The be, br, and rr cuts are shown in the solid lines of colors 
in blue, green, and red, respectively. We ignore the redshift interval of 0.4 < z < 0.45 due to its high sample incompleteness. (A color 
version of this figure is available in the online journal.) 
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Fig. 8. — 1(T scatter (width) of the red sequence galaxies as a 
function of redshift for different magnitude intervals. The solid 
lines are the measured scatter, and the dotted lines arc the intrinsic 
scatter by excluding the photometric errors. (A color version of this 
figure is available in the online journal.) 

readily modeled (Wake et al. 2008; Tojeiro et al. 2012b). 
Comparing to such predictions allows a rough determina- 
tion of the extent of evolution in the red galaxy samples. 

We construct three such samples for the red galaxies, 
which may be expected to resemble a passively evolving 
population, with fixed low, moderate and high number 
densities. For each sample, the fixed number density 
is achieved by finding a (redshift-dependent) luminosity 



threshold AIi{z) and selecting all galaxies with luminos- 
ity above this threshold, as shown in Figure 11. The 
luminosity and redshift ranges are chosen to reduce the 
sample incompleteness caused by the selection cuts. The 
low, moderate, and high number density samples have 
n{z) of 0.4 X 10-4, 1.2 x lO""^, and 2 x IQ-^ h^M-pc'^ , 
respectively. 

As seen in Figure 11, the fixed number density thresh- 
olds correspond globally to rough luminosity thresholds, 
decreasing with increasing number density as expected. 
The luminosity thresholds Mi{z) stay roughly constant 
for the red galaxies at z > 0.48. Since the luminosities 
in our study have been k-l-e corrected, this result implies 
that the stellar population in these red galaxies evolves 
passively. The drop below z < 0.48 is likely caused by 
the incompleteness of galaxies at lower redshift due to 
the CMASS selection cuts, as discussed in Section 2.2. 

The projected 2PCFs for the fixed number density 
samples of red galaxies are presented in Figure 12. The 
top panels show that the clustering strength is inversely 
proportional to the number density, such that galaxies 
in the lower number density samples are more clustered. 
These results are generally consistent with the luminosity 
dependence discussed in Section 3.2, as the different val- 
ues of n{z) effectively act as different luminosity thresh- 
olds. The shapes of the 2PCFs for different n{z) sam- 
ples are similar. The bottom panels compare the 2PCFs 
of the same n(z) samples at different redshift intervals. 
We find that the 2PCFs at different redshifts have simi- 
lar clustering strength, again consistent with the results 
for the luminosity dependence of red galaxies (Figure 4, 
bottom-right panel). We discuss the implications of these 
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Fig. 9. — Measurements of the 3D 2PCF ^{rp,r^) for the blue (left panel) and red (right panel) galaxies in the whole CMASS sample, 
defined using the color cut in Equation 7. Contour levels shown are 5(rp, r-,r) = [0.5, 1, 2, 5, 10, 20]. 



results next. 

3.5. Galaxy Bias 

With the measured galaxy 2PCFs and the theoretical 
matter 2PCF for the cosmology adopted in this paper, we 
can infer the galaxy linear bias factor b from the square 
root of the ratio between the galaxy and dark matter 
2PCFs. The evolution of the linear bias factor can pro- 
vide hints about the evolution of galaxy samples. In 
this subsection, we present the results on galaxy linear 
bias factor 6 as a function of galaxy luminosity, color, 
number density, and redshift, and discuss the possible 
implications on galaxy evolution. Since blue galaxies in 
the CMASS sample are generally far from complete, we 
will focus our discussion on the bias evolution of the red 
galaxies. 

At each redshift, for each galaxy sample, we estimate 
the linear bias factor b by taking the square root of the 
ratio between the measured projected galaxy 2PCF and 
the non-linear dark matter projected 2PCF computed 
at the corresponding redshift, where the latter is cal- 
culated using a modified halof it model (Smith et al. 
2003) with the Eisenstein & Hu (1998) power spectrum 
parametrization. More specifically, we fit the ratio of 
galaxy and matter Wp{rp) with a single parameter b on 
scales 3/i~^Mpc < < 25/i~^Mpc, using the full co- 
variance matrix of the galaxy Wp^Vp). 

Fry (1996) shows that the passive evolution prediction 
for the linear bias factor b(z) follows 

+ (19) 

where D(z) is the linear growth factor at redshift z and 
6o is the bias factor at z = {D{0) = 1). From this 
relation, the redshift evolution of galaxy 2PCF for the 
passively-evolving population can be expressed as 

C(z) = [b{z)D{z)]^,n{0) = [D{z) + {bo-l)]^U{0), (20) 

where Cm(0) is the matter 2PCF at z = 0. Here we 
use the word "passive" to mean that during the evolu- 
tion, each galaxy in the sample keeps its identity and 



there is no merger or disruption that changes the pop- 
ulation. For the CMASS sample considered in this pa- 
per, the galaxy bias factor is usually greater than unity. 
Therefore, according to the above two equations, for 
passively-evolving galaxies, we expect that with decreas- 
ing redshift the amplitude of 2PCF increases while the 
bias factor decreases. 

Figure 13 shows the bias redshift evolution for luminos- 
ity bin samples (left) and for the fixed number density 
samples (right). The bias factors are measured for non- 
overlapping redshift bins. The fitted bias factors sup- 
port our results in the previous subsections that more 
luminous (or lower density) samples are more strongly 
clustered. The dotted lines in Figure 13 are the best-fit 
Fry (1996) relation for each sample, with bo as the single 
fitting parameter. For the two luminosity-bin red galaxy 
samples, we find that they roughly follow the passive evo- 
lution prediction. Strictly speaking, each luminosity bin 
sample does not conserve the number density at different 
redshifts. So by definition it is not a passively-evolving 
population. However, these number density differences 
may be accounted for by slightly changing the luminosity 
thresholds of the luminosity bin sample at each redshift 
(e.g., as a result of an imperfect k-fc correction). The 
measured bias would not be sensitive to such an adjust- 
ment. In such a sense, comparing their bias evolution to 
the Fry relation can still be meaningful. 

For both the luminous and faint samples, there is sug- 
gestive evidence that the clustering at intermediate red- 
shifts (z = 0.575 for the luminous sample, and z = 
0.525 for the faint sample) is slightly weaker than that 
from the best-fit passive evolution. Similar deviations 
from passive evolution were found in other works (e.g.. 
White et al. 2007; Wake et al. 2008; Sawangwit et al. 
2011). However, the large measurement errors in our re- 
sults make the deviation only at about la level, limiting 
our ability for a solid conclusion. 

If this deviation is robustly established, with more ac- 
curate future measurements using CMASS data over a 
larger survey area, it would imply a significant contri- 
bution from processes that break passive evolution, such 
as feedback from active galactic nuclei shutting off star 



14 



Guo ct al. 



10" 



" 1000 



I I 1 1 1 III 1 — I I I IIIJ 



22.5<M,<-21.5 



0.45<z<0.55 



100 r 



10 r 




I I I I I I ' 




-r reddest 
:: redseq 
II green 
blue 



0.1 1 10 100 0.1 1 10 100 

r [Mpc/h] r [Mpc/h] 



10" 



? 1000 r 



o 
a 



100 



10 r 



1 r 



-23<M,<-22 
redseq 



0.5<z<0.6 

0.55<z<0.65 

0.6<z<0.7 



I I I I I I L_I_I_LLLL 



T—rrrmj 1 — i i 1 1 iii| 1 — i i i i iii| 1 — i i 1 1 1 1± 



J X 

i 



-23<M,<-22 
reddest 



0.5<z<0.6 

0.55<z<0.65 

0.6<z<0.7 



I I I I I I 



0.1 1 10 100 0.1 

[Mpc/h] 



10 100 



Tp [Mpc/h] 



Fig. 10. — Color dependence of the projected 2PCF Wp{rp). Top panels display the color dependence for the —22.5 < Mi < —21.5 sample 
at two different redshift intervals. The redshift evolution of uip(rp) for the redseq and reddest subsamples of —23 < Mi < —22 galaxies is 
shown on the bottom. (A color version of this figure is available in the online journal.) 



formation, disruption of satellites in massive halos, and 
mergers of galaxies (Bell et al. 2004; Faber et al. 2007; 
Skelton et al. 2009). The signatm-e may be related to 
the overall migration of blue galaxies to the red sequence 
(Martin et al. 2007), in which case, it would indicate ap- 
preciable migration by z~0.55 — 0.6, consistent with the 
prediction of Faber et al. (2007) . A sophisticated model 
is needed to disentangle the contributions from the dif- 
ferent evolutionary processes. 

The more reliable samples to study the passive evo- 
lution are the ones with fixed number densities, as de- 
scribed in Section 3.4. The results for our three samples 
are shown in the right panel of Figure 13. The low n{z) 
sample appears consistent with passive evolution in the 
redshift range 0.45 < z < 0.65, within the (large) error 
bars on the measurements. This behavior is similar to 
the —23 < Mi < —22 sample. In fact the low n{z) sample 
is close to a luminosity-threshold sample of Mi < —22, 
as shown in Figure 11. For the two samples with higher 



n(z), their bias evolution is consistent with passive evo- 
lution for the smaller redshift ranges probed, which is in 
agreement with the conclusions of Tojeiro et al. (2012b). 

Finally, we present the dependence of the bias factor 
on luminosity in Figure 14 for all galaxies and for the 
red galaxies only. Generally, more luminous and redder 
galaxies at higher redshifts have larger bias factors. For 
the fainter samples, the bias factors are similar in the 
two cases, since at the faint end red galaxies dominate 
the CMASS sample (the majority of the faint blue galax- 
ies are excluded by the selection cuts). The observed 
dependence of galaxy bias factor on galaxy luminosity 
is broadly similar to that for local galaxy samples (e.g., 
Norberg et al. 2001, Zll), but it is non-trivial to compare 
in detail due to the many differences in sample selection, 
redshift, k-l-e corrections and magnitudes. 

We fit the bias-luminosity relation with a commonly- 
used simple functional form (Norberg et al. 2001; 
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Fig. 11. — Construction of fixed number density samples for the red galaxies, as defined by Equation 7, corresponding to low, moderate, 
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distribution of CMASS galaxies. (A color version of this figure is available in the online journal.) 



Zchavi ct al. 2005b), 

b/hp = ci + C2L/Lp 



(21) 



We define Lp as the mean luminosity of galaxies in the 
faintest luminosity bin. This sample of galaxies has 
b = bp, so by construction, ci = 1 — C2. We fit this 
functional form to the luminosity-dependent bias mea- 
surements at 0.43 < z < 0.55, finding C2 = 0.33 for 
all galaxies and C2 = 0.35 for the red galaxies (shown 
as the red dotted curves in Figure 14). At higher red- 
shifts, 0.55 < z < 0.7; the bias factors for all galaxies 
show a decrease compared to the lowcr-redshift relation 
due to the inclusion of more blue galaxies at high red- 
shifts. In contrast, the bias factors for the red galaxies 
globally increase with redshift, as expected from passive 
evolution. The dashed curve in the right panel shows the 
low-redshift relation shifted to high redshift according to 
the Fry relation prediction. 

4. CONCLUSION AND DISCUSSION 

In this paper, we measure the luminosity and color de- 
pendence of the galaxy 2PCFs based on -260, 000 BOSS 
CMASS DR9 galaxies over a —3, 300 deg^ survey area in 
the redshift range of 0.43 < z < 0.7 and study the impli- 
cations on galaxy formation and evolution. 

We first measure the 2PCF for the entire sample. If 
approximated by a power-law, the 2PCF has a correla- 
tion length of To = 7.93 ± 0.06ft-~^Mpc and a slope of 
7 = 1.85 ± 0.01, consistent with the measurements pre- 
sented by White et al. (2011) andNuza et al. (2012). We 
also construct color and luminosity subsamples. To re- 
duce the influence of sample incompleteness caused by 
the target selection criteria of CMASS galaxies, we care- 
fully account for the selection cuts in the color-luminosity 
distribution of galaxies at each redshift interval in defin- 
ing our subsamples. These subsamples in certain color 
and luminosity bins arc close to complete and volume- 
limited. In order to compare the clustering of galaxy pop- 
ulations at different redshifts, we perform k-l-e correction 



to the magnitudes and colors of the galaxies using FSPS 
models. Such corrections are unavoidably dependent on 
the specific stellar evolution models used. However, as 
shown in Appendix A, our general results on the cluster- 
ing analysis would not be significantly affected. 

We find that for all redshift intervals probed, more 
luminous galaxies are more strongly clustered, consis- 
tent with previous studies for galaxies at different red- 
shifts, such as SDSS-I/II main sample galaxies at z — 
0.1 (Zehavi et al. 2005b, 2011), SDSS LRG galaxies at 
z - 0.35 (Zehavi et al. 2005a), and DEEP2 galaxies at 
z ^ 1 (Coil et al. 2006). At each redshift, the large- 
scale galaxy bias factor of CMASS galaxies shows a lin- 
ear dependence on galaxy luminosity, similar to that for 
lower redshift galaxies (Norberg et al. 2002, Zll), but 
with different coefficients in the bias-luminosity relation. 
We divide galaxies globally into a blue and a red popu- 
lation. For each population, we find a similar clustering 
trend - an increasing clustering strength with luminos- 
ity. For blue galaxies, our results are in line with that 
of Zll for SDSS galaxies. For red galaxies, Zll find 
that both the most luminous and faintest galaxies ex- 
hibit stronger small-scale (< 2/i~^Mpc) clustering than 
the samples of intermediate luminosity, which can be ex- 
plained as a large satellite fraction in the faintest sample 
and high mass of host halos for the most luminous sam- 
ple. Because CMASS selects mostly luminous galaxies, 
we are not able to investigate the trend towards faint red 
galaxies, but for the luminous red galaxies, our results 
agree with that in Zll. 

We further investigate the dependence of clustering 
on galaxy color, using finer color cuts. For fixed red- 
shift and luminosity, we find that redder galaxies exhibit 
stronger clustering, similar to the trend found for the 
SDSS main galaxies (Zehavi et al. 2005b, 2011; Li et al. 
2006). Interestingly, such a trend exists even within the 
red sequence, consistent with the finding of Zehavi ct al. 
(2005b, 2011). The trend is different from that of DEEP2 
galaxies in Coil et al. (2008), where no clear color depen- 
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Fig. 12. — Projected 2PCFs, uip(rp), for fixed number density samples of CMASS red galaxies. The top panels show the dependence 
on number density for two redshift bins, and the bottom panels show the redshift dependence for the moderate and low number density 
samples. (A color version of this figure is available in the online journal.) 



dence is seen across the red sequence. If this difference 
is caused by galaxy evolution, it implies that the color 
dependence in the red sequence emerges during the red- 
shift range of 0.7 < z < 1.0. The emergence of the 
dependence may signal the contribution of substantial 
amounts of mergers and inflow of blue galaxies to the 
buildup of the red sequence. 

We also construct subsamples of red galaxies with fixed 
number densities by applying redshift-dependcnt lumi- 
nosity thresholds, and compare their clustering with the 
theoretical prediction of passively-evolving galaxies (Fry 
1996). We find that the evolution of the large-scale 
galaxy bias factors for all the three CMASS subsamples 
considered in this paper are consistent with that from 
the Fry relation, within the relatively large uncertainties 
in the measured bias factors, which suggests that the 
red galaxies in the CMASS sample roughly follow pas- 
sive evolution from z — 0.7 to 0.45. In contrast, from 
HOD modeling of clustering of red galaxies in NDWFS 



(Brown et al. 2008), White et al. (2007) found that pas- 
sive evolution from z ~ 0.9 to z ~ 0.5 would predict too 
many satellite galaxies in high-mass halos and concluded 
that about one-third of these satellites must have expe- 
rienced merging or disruption. The apparent discrep- 
ancy between our result and that in White et al. (2007) 
can be explained by the difference in the number densi- 
ties of galaxy samples. The NDWFS samples analyzed 
in White et al. (2007) have a constant comoving num- 
ber density of n{z) ~ 10~'^/i'^Mpc~'^, about one order 
of magnitude higher than the ones we study here. Thus, 
their constant number density samples include fainter red 
galaxies (which have a larger contribution from satel- 
lite galaxies). In contrast, galaxies in our more lumi- 
nous samples are predominantly luminous central galax- 
ies that roughly follow passive evolution (but see also 
Wake et al. 2008; Sawangwit et al. 2011). These results 
seem to support a scenario in which mergers and disrup- 
tion play an important role for the evolution of low-mass 
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Fig. 13. — Linear bias factor, b, as a function of redshift for the CMASS red galaxies in different luminosity samples (left panel) and different 
fixed n{z) samples (right). The symbols represent the measured galaxy bias b{z) by fitting Wp{rp) over 3 h~^Mpc < rp < 25 h^'^Mpc relative 
to the theoretically predicted one for dark matter. The dotted lines are the best-fit passive evolution predictions by fitting b(z) using the 
Fry (1996) relation. The goodness of fit x^/^o/ is also given for each set of samples. 



3.5 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 
All galaxies 

O 0.43<z<0.55 

□ 0.55<z<0.7 



n — I — I — I — I — I — I — I — I — I — I — I — I — I — I / 1 . I 

/ .•■ 
/ .' 



3 - 



2.5 - 



2 - 



••t] 



.35' 



_i I I i_ 



_1 I I L_ 



_I I I L_ 



Red galaxies 
O 0.43<z<0.55 
□ 0.55<z<0.7 



[] / 
/ /' 

/ 

/ .■■ 



/ 



/ 



/ . 



_l I I L_ 



_l I I L_ 



_l I I L_ 



-21.5 -22 -22.5 -23 -21.5 -22 -22.5 -23 



Fig. 14. — Linear bias factor, b, as a function of luminosity in two redshift ranges for all CMASS galaxies (left) and for the red galaxies 
only (right). The dotted curves arc the bias-luminosity relation. Equation 21, with C2 = 0.33 for all galaxies and C2 = 0.35 for the red 
galaxies, for the lower redshift range. The dashed curve in the right panel is the low-rcdshift relation shifted to the higher-redshift range 
according to the passive evolution prediction. 



red galaxies. 

Our investigation of the color-luminosity distribution 
at each redshift reveals two notable trends in the width 
of the red sequence. The red sequence becomes narrower 
towards the high-luminosity end, and it becomes nar- 
rower towards lower redshifts. Similar results are also 
seen in galaxies at both lower and higher redshifts (e.g., 
BeU et al. 2004; Skibba & Sheth 2009; Whitaker et al. 
2010). The color scatter in the red sequence reflects 
the distribution of the ages of stellar population, dust 



extinction, and metallicity. At a given redshift, fainter 
galaxies show a more diverse distribution of these quan- 
tities, leading to a wider distribution in color. Passive 
evolution makes galaxies redder and largely reduces the 
color difference caused by the distribution of the ages 
of stellar population, leading to a narrower red sequence 
towards lower redshifts. 

The inferences in this paper about the evolution of 
CMASS galaxies from the measured color and luminos- 
ity dependent clustering are still based on simple clus- 
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Fig. 15. — Color-magnitude diagram for CMASS galaxies using the FSPS and Mil stellar evolution models, for two rodshift intervals of 
0.5 < z < 0.55 and 0.6 < z < 0.65. The red contour lines are for the FSPS model and the blue ones are for the Mil model. The green line 
is our adopted color cut using the FSPS models. The points arc the galaxies from the CMASS Sparse Sample (see text). (A color version 
of this figure is available in the online journal.) 

tering models and interpretations. A further, natural 
step to interpret these results is to perform HOD mod- 
eling of our measurements, which will allow us to better 
study galaxy formation and evolution by incorporating 
knowledge about the dark matter halo formation and 
evolution. We expect that improved measurements from 
larger BOSS samples in the future and detailed HOD 
modeling will greatly advance our understanding of the 
evolution of massive galaxies. 
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APPENDIX 

A. DIFFERENT STELLAR EVOLUTION MODELS 

In this paper, we correct for the k-f-e effects using the FSPS model, as mentioned in 2.2. Tojeiro et al. (2012b) 
compare the FSPS model and the stellar evolution model of Maraston & Stromback (2011; Mil), and conclude that 
both models provide similar star formation histories with similar mass-weighted ages. Tojeiro et al. (2012b) also 
compare the large-scale clustering using both models, and find that they produce consistent results. We show in 
Figure 15 the color-magnitude diagram obtained using the two different models at two redshift ranges. Contour lines 
for the FSPS model are shown in red and for the Mil model in blue. The green line is our proposed color cut (using 
the FSPS model). The two models predict similar k-l-e corrections for z < 0.55, but at higher redshifts the Mil 
model appears to produce more luminous and bluer galaxies. The overall shapes of the distributions for the two 
models, however, are quite similar. In particular, using alternatively the Mil model would not tilt the slope of the red 
sequence. Using the Mil model, the detailed color cuts would be changed accordingly but the clustering dependence 
on the luminosity and color are not expected to change. 

From the CMD, it appears that CMASS galaxies have a tail of faint extremely-red galaxies, which is different from 
the properties of the SDSS main sample (see, e.g., Zll). This appearance, however, is partly caused by the selection 
cuts, which remove most of the faint blue galaxies. To clarify this effect, we show in Figure 15 the galaxies from the 
BOSS CMASS Sparse Sample, which includes fainter and bluer galaxies by extending the sliding cut (Equation 3) to 

tcmod < 20.14 + 1.6(di - 0.8) (Al) 

(Padmanabhan et al. 2013). These faint blue galaxies are sparsely sampled to yield approximately five objects per 
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Fig. 16. — Diagonal jackknifc errors on the projeeted 2PCF, Wp, estimated from moelc catalogs (left-hand side) and the real data (right), 
using different number of jackknife samples. In the left panel, the solid curve shows the fractional errors estimated from the variance among 
Wp measurements in 100 mock catalogs (Manera et al. 2012). The symbols are the average (over the 100 mocks) of the fractional jackknife 
errors in each catalog. The errorbars plotted reflect the variation among the jackknife estimates across the 100 mocks. The number of 
jackknifc samples in each mock ranges from 50 to 200, as indicated by the color of the symbols. On the right we show the fractional 
jackknife errors estimated from our "one realization" of the actual CMASS data, for different number of jackknife samples. (A color version 
of this figure is available in the online journal.) 

square degree. It is clear that the red sequence has more triangular distribution in CMASS, as a result of the increase 
in the width of the color distribution for fainter galaxies (Figure 8 and discussion thereof). A similar shape of the 
red sequence is also observed at higher redshifts (see e.g., Figure 1 of Whitaker et al. 2010). The narrow red sequence 
shown in Zll is caused by the fact that their CMD is for the flux-limited main sample galaxies. Many fainter galaxies 
at higher redshifts are excluded by the faint flux limit. If shown for a volume-limited sample, the "triangle" shape of 
the red sequence is more apparent (see e.g., Figure 2 of Skibba & Sheth 2009). Since this Sparse Sample does not have 
the same selection as other CMASS galaxies, we do not include it in our analysis, but these galaxies can be useful for 
studying evolution of blue galaxies that have smaller stellar masses. 

B. JACKKNIFE ERROR ESTIMATES 

We use the jackknifc resampling method to construct the error covariance matrices (Equation 12). In principle, it 
is preferable to derive covariance matrices from large numbers of realistic mock catalogs, each matching the observed 
galaxy properties, survey geometry, and selection functions. However, our tests below demonstrate that the jackknife 
error estimates perform quite well and arc sufficient for our purposes. Moreover, it is a far more practical tool 
when working with many subsamples of different clustering properties. We do not currently have available mocks 
with suitable modeling of the galaxy luminosity function and with the correct galaxy distribution on small scales. 
The large set of mock catalogs used by Anderson et al. (2012) are constructed by populating dark matter halos 
in simulations with galaxies according to the HOD model fitted to the redshift-space 2PCF ^(s) on large scales 
{30 h~^M-pc < s < 80/i~^Mpc). The measured large-scale 2PCFs of galaxies are reasonably reproduced in these 
mocks, while the small-scale 2PCFs are not matched (see Manera et al. 2012, for details). Although the small scale 
2PCFs measured from these mocks generally deviates from the real data, the mocks can still be used to evaluate the 
validity of the jackknife method, and the appropriate number of jackknife subsamples to use. 

In the left panel of Figure 16, we compare the fractional diagonal errors of Wp from variations among 100 mock 
catalogs of Manera et al. (2012) (solid curve) and those from applying the jackknife resampling method to the mocks 
with different number of jackknife samples (symbols with different colors). The jackknife error estimates show excellent 
consistency with those estimated from variations among mocks, especially for scales less than a few /i~^Mpc. Around 
10ft,~^Mpc, the jackknife method slightly overestimates the errors (by ~15%). The jackknife errors also appear to be 
insensitive to the total number of jackknife samples used, with a general convergence for 100 or more jackknife samples. 
We have also checked the off-diagonal elements in the covariance matrices and find that the jackknife estimates are 
generally consistent with the variations-among-mocks estimates, although there are somewhat larger fluctuations of the 
jackknifc estimates from mock to mock. All the results are consistent with those found and discussed in Zehavi et al. 
(2005b). 

Norberg et al. (2009) have performed a comprehensive study of a variety of error estimators for dark matter corre- 
lation functions in A'^-body simulations, comparing estimates such as jackknife and bootstrap to those derived from 
multiple independent catalogs. They find good agreement between jackknife and external estimates for the variance in 
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Fig. 17. — The covariance matrix of Wp for the whole CMASS sample, estimated with 100 jackknife samples, normalized by the diagonal 
elements of the matrix. The scale on the right shows the color scheme, representing the level of covariance on the different scales. (A color 
version of this figure is available in the online journal.) 

Wp(rp) on large scales, consistent with our results in Figure 16. However, Norberg et al. (2009) find that the jackknife 
method significantly overestimates the errors in Wp on small scales (e.g., by 40% at rp < l/i~^Mpc), while we do 
not see such a large difference in our results. It is worth noting that we differ in our implementation of the jackknife 
method. While Norberg et al. (2009) estimate the jackknife errors by dividing the simulation boxes into N subvolumes, 
we divide the sample into N jackknife subsamples of equal area with the same radial selection. Also, the mocks used 
in Norberg et al. (2009) have smaller volume than the ones used in this paper, which might lead to larger uncertainties 
in their error estimates. It is possible that the differences in applying the jackknife method and in the uncertainties 
can explain the apparent difference between our results and those of Norberg et al. (2009). 

The right panel of Figure 16 shows the jackknife error estimates for the entire CMASS sample. It is encouraging 
that the overall shape and magnitude of the fractional errors from the real data are in good general agreement with 
those from the mocks, even though the mocks do not intend to match the small-scale clustering. Since we do not 
see significant variations in the error estimates with 100 or more jackknife samples, we choose N = 100 jackknife 
samples for the error estimation in this paper. With an effective area of about 3, 300 deg^ for the CMASS DR9 sample, 
each jackknife sample then has an area of ~33deg^, corresponding to about 2.1 x lO'' (/i~^Mpc)^, large enough for 
measuring the 2PCF in the range presented in this paper (< 50 /i~^Mpc). The normalized jackknife covariance matrix 
with N = 100 samples, for Wp of the whole CMASS sample, is shown in Figure 17. While the correlation coefficients 
for the off-diagonal elements on small scales (< 2 h^^Mpc) are mostly below 0.3, the errors of Wp on large scales are 
highly correlated. Therefore, when fitting Wp for scales larger than 2 h^^Mpc, the full covariance matrix, not just the 
diagonal elements, should be taken into account. 
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